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Topological inhomogeneity gives rise to spectral anomalies that can induce Bose-Einstein Con- 
densation (BEC) in low dimensional systems. These anomalies consist in energy regions composed 
of an infinite number of states with vanishing weight in the thermodynamic limit (hidden states). 
Here we present a rigorous result giving the most general conditions for BEC on complex networks. 
We prove that the presence of hidden states in the lowest region of the spectrum is the necessary 
and sufficient condition for condensation in low dimension (spectral dimension d < 2), while it is 
shown that BEC always occurs for d > 2. 



The recent experimental evidence for Bose-Einstein condensation in real systems [1] has stimulated an increasing 
amount of theoretical works to find the most general conditions inducing such phenomenon in new experimental 
setups. A key problem in this analysis is the influence of geometry on the physical behaviour of bosonic models. In 
the classical framework of the ideal Bose gas this influence is encoded in the dependence on the space dimension d of 
the system: indeed BEC occurs if and only if d > 3 [2] . 

Recently, the interest in quantum devices such as Josephson junctions, together with the possibility of combining 
them in complex geometrical arrangements, has stimulated the study of bosonic models on general discrete structures 
[3]. There, the problem of defining an effective dimension describing large scale topology has been successfully solved 
by the introduction of spectral dimension d [4,5], which can be experimentally measured [6] and rigorously defined by 
graph theory [7]. On the other hand, one expects the influence of topology to be richer and more complex on discrete 
structures, due to the possible relevance of local geometrical details, in addition to the large scale structure described 
by dimensionality. 

In this direction, recent works on BEC on inhomogeneous networks [8] have put into evidence that strong inhomo- 
geneities can give rise to condensation at finite temperature even in a low-dimensional system as a comb lattice (d = 1, 
see figure 1). This phenomenon arises from a peculiarity of the spectra of inhomogeneous networks we shall refer to 
as hidden spectrum, consisting of an energy region filled by a finite of infinite number of eigenvalues which do not 
contribute to the normalized spectral density in the thermodynamic limit. Hidden spectra don't usually affect bulk 
thermodynamic quantities but, as we shall show in the following, can have dramatic effect when bosonic statistics 
forces the macroscopic occupation of a single quantum state. 

Here we give a general mathematical definition of hidden spectra and we prove that the necessary and sufficient 
condition for condensation when d < 2 is the presence of hidden states in the lowest region of the energy spectrum, 
while BEC at finite temperature always occurs for higher dimensional systems (i.e. when d > 2). 

A generic discrete network is naturally described by a graph which is a countable set V of vertices (sites) i connected 
pairwise by a set E of unoriented edges (links) (i,j) = (J,i)- We will call nearest neighbours two vertices joined by 
an edge. The graph topology is algebraically described by its adjacency Ay = 1 if is a link of the graph and 
Aij = otherwise. The coordination number z% = Y]j is the number of nearest neighbours of the site i. A real 
geometrical structure always has a finite maximum number of nearest neighbours at any site. To take into account 
this condition in the thermodynamic limit we introduce the uniform boundedness condition: maxjZj < oo on the 
coordination numbers. A path in G is a sequence of consecutive links {(i, k), (k, h), (n, m), (ro, k)} and a graph 
is said to be connected if for any two point there is always a path joining them. In the following we will consider 
only connected graphs; disconnected graphs can be reduced to non-interacting connected components which can be 
studied separately. Every connected graph is endowed with an intrinsic metric generated by the chemical distance 
7"jj, which is defined as the number of links in the shortest path connecting sites i and j. 

The most general Hamiltonian for non interacting particles on G is: 
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H = J2 h H a \ a 3 (!) 
i,jev 

where a\ and dj are the creation and annihilation operator at site i. The bosonic nature of the particles is introduced 
through the usual commutation relations [dj, at] = Sij. The Hamiltonian matrix /i^ is defined by: 

hij = tij + S i: jVi (2) 

The term tij describes hopping between nearest neighbours sites and it is directly related to the topology of the graph. 
Indeed tij ^ if and only if G E i.e. Aij = 1. The diagonal term Vt takes into account a potential at site i and 
both terms must satisfy a a boundedness condition: < k < \Uj\ < K and < c < \Vi\ < C. 

To introduce the thermodynamic limit and analyze the conditions for Bosc-Einstein condensation on the graph, we 
begin by studying models restricted to the Van Hove sphere S ry0 of center o and radius r. This is defined as the set 
of the vertices whose distance from o is equal or smaller than r. We will call N rj0 the number of site in S _ r . The 
Hamiltonian restricted to the sphere is: 

RSr '° = E h iP°4aj (3) 
i,jev 

s s 
where h^'" — hij if i and j belong to the sphere and h^'° = otherwise. For graphs with polynomial growth, i.e. 

when N r ^ ~ r p for r — > oo, it is possible to show that the thermodynamic limit is independent from the choice of 

the center of the sphere o [9]. In the following we will consider only graphs with polynomial growth, since this is the 

necessary condition to ensure that the structure can be embedded in a finite dimensional space. We will then write: 

H r = H s r , % h r, = h Sr.o and Nr = Nr o _ 

For each finite sphere of radius r the matrix h\j defines a normalized density of states p r (E) which is the sum of 
N r (5-functions S(E — E'Q, where E T k are the eigenvalues of h\j. The density p r (E) will be normalized to 1/N r . 
In the thermodynamic limit, we define p(E) to be the spectral density of the eigenvalues of hij if: 

lim / \p r (E) - p(E)\dE = (4) 

r — >oo / 

Let us define E m = Inf(Supp(p(S))) where Supp(p(£?)) is the support of the distribution p(E). The asymptotic 
behavior of the thermodynamic spectral density in this region is described by the spectral dimension d [5]: 

p(E) ~ {E-E^i- 1 for E^E m (5) 

A hidden region of the spectrum is an energy interval [^1,^2] such that [Ei,E 2 ] n Supp(p(E)) = and 
linv^oo N^ Ei E ^ > 0, where N£ Ei E ^ is the number of eigenvalues of h\j in the interval [Ei,E 2 ]- Notice that in 
general N^ Ei E ^ can diverge for r —> 00 and the eigenvalues can become dense in [E\, E 2 ] in the thermodynamic limit. 
Therefore this condition not only includes the trivial case of discrete spectrum but is far more general; an interesting 
example of this new kind of behaviour has been observed in the comb lattice [8] (figure 1), where the hidden region 
of the spectrum is filled by an infinite number of states. 

We now define the lowest energy level for the sequence of densities p r {E), setting Eq = Inffc(.E£) and Eq = 
linv^oo Eq. In general, E < E m . If E < E m , then [E , E m ] is a hidden region of the spectrum, which will be called 
hidden low energy spectrum. 

In the following we will consider models at fixed fillings / = N/N r , where TV represents the number of particles in 
the system. In the macro-canonical ensemble the equation that determines the fugacity z in the thermodynamic limit 
is: 

/ = lim f (6) 

Setting Eq = we have that < z < 1. 

The integral in the equation (6) can be divided into two sums, the first ones considering the energies smaller than 
an arbitrary constant e and the second the energies larger than e: 

f P r (E)dE 1 f p r(E)dE 

J z-^ePE - l Z-> z -i e fJEi _ x + J E>e z -i ef 3E _ 1 V) 
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We define: 



E k <e ^ 

< S E z -l e ,Er _ 1 (8) 

k=0 

as the fraction of particles with energy smaller than e. Bose-Einstein condensation occurs in this systems if it 
exists a critical temperature T c > such that, for any T < T c , n e = limbec n r t > k > 0, for all e > 0; i.e. 
n = lim e _ >0 n e — k > 0. 

From the definition (8), no can be strictly positive only if lim,.-^ z(r) = 1. Indeed, if this limit is smaller than 1 
it follows that: n e < (1 — z)^ 1 linv^oo Nl/N r where is the number of state with energy smaller than e. If the 
ground state is not infinitely degenerate, lim e ^ lim r ^ 00 (iVT/iV r ) = and then n = 0. 

Taking first the limit r — > oo and then e — > in equation (7) we obtain: 

e^o™/^ z^e^ - 1 e^o™J fi>£ z^e^ - 1 

Now, from the boundedness of (z~ 1 e /3 - E — for £7 > e and from the definition (4), the first of the two limits in 
the right hand side of (9) vanishes: 

p{E)dE f p(E)dE 



f = n ° + 1™ L z-we n = n ° + / z -^e i (io) 



where, again, no can be different from only if z = 1. 

Now the integral in equation (10) is an increasing continuous function of z with < z < 1. If the limit: 

Aw-ja/r^r (11> 

is equal to oo we have z < 1 and n = 0. If the limit is finite, / c (/3) is a decreasing function of j3 with lim^oo / c (/3) = 
and lim/3^0 /c(/3) = oo. Then, for a suitable 0c, fc(0c) = /• For > C , (i.e. T < T c ) z = 1 and n = f - f c (0) > 
while for < C , (i.e. T > T c ) z < 1 and n = 0. 

From the divergence or finiteness of the limit (11) one obtains the most general conditions for the occurrence of 
Bose-Einstein condensation. 

First, if = E < E m (i.e. the system present a low energy hidden spectrum) the limit (11) is finite and there is 
Bose-Einstein condensation at finite temperature: 

f°W * ^m—i J p^ dE - (12) 

On the other hand, when E = E m the value of the limit (11) is determined by the spectral dimension d. Indeed if 
d > 2: 

where 5 and c\ are suitable constants. Therefore in this case Bose-Einstein condensation occurs at finite temperature. 
For d < 2 we have: 

s z Cl Ei- l dE 



M>**J o W T^ = ° (14) 

and there is no Bose-Einstein condensation. When d = 2 we have to consider the logarithmic correction to (5) and it 
is possible to show that the limit (11) diverges. 

This result applies to many different situations. The simplest example is the discretization of the usual Schrodingcr 
equation for free particles on a graph. In this case the Hamiltonian is given by: hij — j^Lij. Lij is the Laplacian 
operator on the graph = ZiSij — Aij . It can be shown that Eq = E m and therefore the occurrence of BEC is 
determined by the behavior of the spectrum of at low eigenvalues, described by the spectral dimension as in (5). 
d is known for a wide class of structures [7,10,11]: for lattices it coincides with the usual Euclidean dimension and for 
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exactly decimable graph (the Sierpinski gasket [10] and the T-fractal are illustrated in figure 2) it can be proven [7] 
that d < 2. 

A more important model, relevant for real condensed matter structures, is a pure hopping of non interacting bosons 
on graphs. This has been considered in [8] for the description of the Joscphson junction arrays in the weak coupling 
limit. In this case the Hamiltonian matrix is given by: hij = —tA^. A relevant point about the application of our 
result to real systems concerns the effects of the introduction of a fluctuating local potential. Indeed it can be shown 
[12] that in presence of a hidden spectrum giving rise to BEC, the existence of a condensate is not affected by the 
introduction of a small enough potential. 

Let us now focus on the pure hopping model described by the Hamiltonian = —tAij where the only effects are 
due to topology, to show how the general theorem can be applied to a wide class of discrete structures. Condensation 
at finite temperature due to the presence of low energy hidden states is typical of bundle structures [11]. These graphs 
are obtained by a "fibering" procedure, i.e. attaching the origin site of copy of a graph we will call the "fiber" graph, 
to every point of another graph called "base" . An example is the brush graph, shown in figure 3, where the base (the 
2d lattice) is fibered by a linear chain. On bundled structures the spectrum p(E) is simply given by the spectrum 
of the pure hopping model on the fiber but the inhomogeneities due to the base give rise to a low energy hidden 
spectrum. Moreover, the wave function of the condensate is localized along the base and presents a fast decay along 
the fibers. This result can be in general obtained by first diagonalizing the Hamiltonian of a pure hopping model 
defined on the base and then solving the eigenvalues problem of the Hamiltonian of the fiber with a suitable impurity 
in the origin, whose form depends on the result of the previous diagonalization. (see [8] for a detailed application 
of this techniques to the comb graph). In the case of the brush graph E m = -\/(20) — 2 « 2.47 > E = (in all 
figures the energy is measured in t units and the energy zero has been chosen so that E — 0), so that we can apply 
the general result for low energy hidden spectra. In the figure the continuous line is p(E) and it coincides with the 
density of states of the hopping model on a linear chain (the fiber). The dotted line represents the hidden spectral 
region, which is filled continuously by the hidden states. Since on the comb graph the fraction of states in the hidden 

1 /3 

spectral region goes to zero as 1/N r for r — > oo, in order to obtain the dotted lines we have to normalize p r (E) 

2/3 

dividing by N r and not by the total number of sites in the sphere. 

Graphs with constant coordination number Zi are typical examples in which the pure hopping model do not present 
hidden low energy regions (Eq = E m ). This is due to the fact that for this class of graphs the spectrum of the model 
can be obtained from that of the Laplacian matrix by a shift of the zero in the energy. The existence of BEC is 
therefore determined by the spectral dimension d of the graph. This parameter can be exactly calculated for a wide 
class of discrete structures. On lattices, d — d is the usual Euclidean dimension and one recovers the classical result 
for BEC on translation invariant structures. The d-dimensional Sierpinski gaskets [10] (see figure 2) are example of 
exactly decimable graphs with constant coordination number. Since for these graphs d < 2, our general result proves 
that in these cases no Bose-Einstein condensation occurs. 

A fundamental property of the low eigenvalues spectral density of the Laplacian matrix is its independence from 
the local details of the graph, i.e. gemoetrical universality [5]. The value of d is not changed under a wide class of 
transformations, called isospectralities, which can strongly modify the geometry of the graph. A simple consequence 
of this property is that if we consider the pure hopping model on a graph with constant coordination number, which 
differs from a graph of known dimension d by an isospectrality, BEC occurs only if d > 2. An example of this 
behaviour is given by the ladder graph (figure 4) which can be obtained from the linear chain by the addition of 
finite-range links, which is one of the simplest cases of isospectralities On this structure then there is no low energy 
hidden spectrum, d = 1 and the pure hopping model does not exhibit BEC. 

An important property of d is that the dimension of the graph obtained as a direct product of two graphs is the sum 
of the dimensions of the original structures. An example is illustrated in figure 5, where we show the direct product 
of a linear chain and a well known fractal, a Sierpinski gasket. In this case the coordination number Zi is constant. 
Applying this property of bard it is easy to show that d = 1 + 21n(3)/ ln(4) > 2, and therefore one immediately infers 
that BEC occurs on this graph. Interestingly, in this case we do not know the analytical form of the spectrum of the 
Hamiltonian (the spectral density in the figure is obtained by a numerical diagonalization), nevertheless we can prove 
the existence of BEC from the general theorem and from the properties of d. 
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FIG. 1. The comb graph and the spectrum of the pure hopping model (hij = —tAij) defined on this structure. The hidden 
spectrum is represented with dotted lines. In this case the spectrum is obtained from an exact calculation in the thermodynamic 
limit. The energy is measured in units of t and the zero has been chosen so that Eo = 0; the density is plotted in arbitrary 
units (the hidden spectrum and p(E) can be normalized to 1 dividing respectively to the number of states belonging to each 
spectral region). 




FIG. 2. Two classical examples of exactly decimable graph: the 3d Sierpinski gasket (d = 21n(4)/ln(6) < 2) and the t-fractal 
(d = 21n(3)/ln(6) < 2). 




FIG. 3. The brush graph and the spectrum of the pure hopping model on this graph (obtained from the diagonalization of 
a structure of 200 x 200 x 200 sites). The dotted lines represent the hidden spectrum which gives rise to condensation even if 
the spectral dimension is 1. 
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FIG. 4. The ladder graph and the spectrum of the pure hopping model obtained from an exact diagonalization in the 
thermodynamic limit. Here there is not hidden spectrum and p(E) — > oo when E — > E m = Eo (d = 1), then there is not 
condensation at finite temperature. 
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FIG. 5. The graph obtained as a product of the Sierpinski gasket and a linear chain. Here the spectrum of the pure hopping 
model is obtained from the diagonalization of a graph of 840192 sites. Here condensation occurs since d = 1 + 21n(3)/ ln(4) > 2 
(indeed the spectral vanishes for E — > E m = Eo). 
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